Spectral Broadening of Radiation from Relativistic Collapsing Objects 
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We discuss light curves and the spectral broadening of the radiation emitted during the finite 
interval of time from the surface of a spherically symmetric collapsing object. We study a simplified 
model of monochromatic radiations. We discuss how one can obtain information about the physical 
parameters of the collapsing body, such as its mass and radius, from the light curves and spectral 
broadenings. 
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Light propagating in the vicinity of astrophysical com- 
pact objects, like neutron stars and black holes, is af- 
fected by the gravitational field. It has been demon- 
strated that the general relativistic effects might be im- 
portant for understanding the features of the radiation 
coming from the neutron star like objects p], 0, Q. The 
gravitational redshift and bending of light rays emitted 
by a compact object affect the form and spectrum of the 
observed signals. For the emission of light from the vicin- 
ity of a black hole these effects are more profound than 
for neutron stars. For example, the line broadening of X- 
ray observed by the ASCA satellite can be explained by 
the strong gravitational effect on the light emitted from 
the accretion disk located near to the central black hole 

Astrophysical black holes are believed to be formed 
as a result of the gravitational collapse of massive stars 
@ . Then it is natural to expect that the radiation emit- 
ted during the gravitational collapse is also affected by 
the strong gravity. For a spherical collapse and contin- 
uous emission of light this effect was studied in details 
0,H@|. Recently the light curves for collapsing objects 
were studied in a slightly different set up [10] assuming 
that the radiation has a profile of a sharp in time pulse. 
Such radiation may occur during the collapse of a star 
when its matter density becomes much higher than the 
nuclear density. Under these conditions hadronic phase 
transitions are expected which may result in sharp- 
in-time emission of massless particles (photons and neu- 
trino) 0. 

In this work, we consider the radiation emitted by a 
collapsing star during a finite time interval and calculate 
light curves and the spectrum of this radiation as seen by 
a distant observer. As in the previous work [10], we adopt 
a simplified model of a freely falling spherical surface and 
assume that the radiation is originally monochromatic. 
But instead of instant radiation, we focus on the radia- 



tion emitted during the finite interval of time. The main 
goal of this study is to analyze how one can extract in- 
formation about the characteristics of a collapsing object 
(its mass and radius) from the observed spectra and light 
curves. 

In the adopted simplified approach, we consider 
a spherical collapse and assume that the points of 
the collapsing surface follows radial geodesies in the 
Schwarzschild geometry 0, OH EH, 03 • Thus the pho- 
tons emitted from the surface propagate to the observer 
at infinity in the background of the Schwarzschild metric 



ds 2 = -fdt 2 



f^dr 2 +r 2 dtt 2 . 



(1) 



Here / = f(r) = 1 — 2M/r, M is the mass of the collaps- 
ing object, and dVt 2 = dO 2 -fsin 2 d(j) 2 is a line element on 
a unit sphere. We use the natural units, c = G = H = 1 
throughout this paper. 

Denote by r the proper time as measured by an ob- 
server comoving with the collapsing surface. We suppose 
that the collapse starts at r = with the initial surface 
radius Ro. We denote by t^ e \r) the Schwarzschild time 
t corresponding to r and take t^ e \r = 0) = 0. In the 
coordinates (t, r, 0, <p) the four- velocity of the collapsing 
surface is ^ = (dt^ /dr,dR/dr, 0,0). We define vj as 
the invariant radial velocity which measures the proper 
length change as measured by the proper time of the ob- 
server at rest at a given radius. One has vi = f~ 1 (R)^ , 
\vj | < 1. For a freely falling surface with the initial radius 
Ro the invariant velocity as a function of R is 



vi 



2M yjl - R/R 



R y/1 - 2M/R 



(2) 



Because of the spherical symmetry the trajectories of 
particles and light are plane, so that we can always put 
= 0. Denote by p^ = (p t ,p r ,p e ,0) the 4-momentum 



2 



of a photon, then E = —p t (the energy at infinity) and 
L = pe (the angular momentum) are constants of mo- 
tion. In the adopted units, where h = 1, E coincides 
with the photon frequency as measured by a far distant 
observer, and p coincides with its wavevector. Instead 
of the angular momentum L we shall use the impact pa- 
rameter defined by I = L/E . The radial momentum p r 
is given by p r = crEZ, where Z(l,r) = y/l — l 2 f(r)/r 2 
. Here and later a denotes a sign function which takes 
the values + and — for a forward (p r > 0) and backward 
(p r < 0) motion of the photon, respectively. 

For a photon emitted at the radius R and propagating 
to the infinity there exists an upper limit for the impact 
parameter l m ax — R/ V ' f{R) determined by the condi- 
tion Z^jnax^R) = 0. Let us consider the emission angle, 
/3, of the radiation as measured by an observer comov- 
ing with the surface. The emission angle is zero for the 
emission in vertical outward direction and ir/2 for the 
tangentially emitted radiation. Since the matter inside 
the collapsing surface is supposed to be opaque, it is nat- 
ural to consider only the light emitted with the emission 
angle, (3 < tt/2. It gives a lower bound on the impact 
parameter for the backw ard emission, It < I < lmax, 
where l T = R/ ^/f(Ro) is determined by the condi- 
tion for tangentially emitted radiation, Z(It,R) = —vi • 
Hence the possible ranges of the impact parameter are 
< I < lmax and It < I < lmax for a forward and 
backward emission, respectively. In this work we con- 
sider only R > 3V3\/l -2M/R M. A discussion of the 
allowed ranges of the impact parameter for the smaller 
radius up to R ~ 2M, can be found in @, Q- 

We choose the direction of the axis 6 = such that a 
photon emitted with (3 = at = propagates along 
the radial direction to a distant observer. For a photon 
emitted from a collapsing surface at the angle to reach 
the distant observer, coincides with the bending angle. 
For a null ray from the collapsing surface when its radius 
is R, the bending angle for a forward-emission, 6+ is 



+ = e(z,i?) = z 



dr 



R 



r 2 Z(l,r) 



(3) 



For a backward-emission, a photon, before it reaches 
the infinity, should first pass through a turning point 
r t (< R). The turning point is determined by Z(Z, r t ) = 0. 
Then we get the bending angle for backward emission 
as 9-(l,R) = 20(Z,r*t) — 0(Z,i?). One can easily check 
that for a tangentially emitted backward photon the an- 
gle 6-(It,R) is greater than tt/2. This means that a 
distant observer can see a part of the 'opposite side' of 
the spherical surface as a result of strong gravity effect. 

Let pft and pft be 4-momentum of a photon emit- 
ted by a collapsing surface and of a photon at infinity 
respectively. Then = —p^v^ is the energy (fre- 
quency) of the photon as measured by a comoving ob- 
server. For the observer at rest at with the 4- velocity, 
v (o)n _ ^ observed energy (frequency) is deter- 
mined by = —p^v(°^. For a given ray with the 



impact parameter Z, which is emitted from the freely- 
falling surface, Eq. ([2]), when its radius is R, the redshift 
factor $ defined as the ratio of the emitted frequency to 



the observed frequency at infinity, $ 

by 



Ae) 



/ v (o) 



1 - aviZ(l,R) 



is given 



(4) 



Consider a light ray with the impact parameter Z emit- 
ted from the collapsing surface at the moment of the 
proper time r, when it has the radius i?(r), Rf(jf) < 
R(r) < Ri(ji), and let t be the time when it reaches the 
distant observer. The first ray reaching the distant ob- 
server is the ray with Z = emitted from a point 'a' at 
Ti (see Fig. 1). We use it as a reference ray. We charac- 
terize the arrival time of other rays by their time delay 
At with respect to the arrival time of the reference ray. 
In Fig. 1 we show the sets of points of emission {i?, 0}, 
which have the same At, by dotted lines. For the forward 
ray the time delay At is given by the following expression 



AM*; T, Ti) = tW(T) - W (Ti) + T(l, R(T)) 

R(t) R *~ 2M 



+Ri 
T(l,R) 



-2Mln 

00 dr ' 



R 



f(r) 



R(t) 
1 



[Z(l,r) 



2M 
- 1 



(5) 
(6) 



Similarly for the backward ray one has 



At-(l;T,n) = &\t) - &\n) + 2T(l,r t ) - T(1,R(t)) 



+Ri + R(t) - 2r t + 2M In 



(R(t) - 2M){Rj - 2M) 
(r t - 2M) 2 



•(7) 



The integrals for and T, Eqs. ([3]) and (|6]) respec- 
tively, can be expressed in terms of the elliptic functions. 
However, for practical calculations it is very convenient 
to use the analytical approximations for these quantities 
in terms of simple elementary functions [10, LL5[ . 

For the monochromatic emission with frequency = 
^, the flux as measured by a distant observer at ro can 
be written as an integral over the proper time, r: 



Fl°\t) 



2?r f 



drl 



(8) 



where J^(Z,r) is the intensity of the radiation (see e.g. 
[10[). In the present paper we focus on the finite-duration 
monochromatic emission, and take the intensity in the 
form li e \l,r) = f(r)X^(l). For simplicity we assume 
that the profile f(r) is a step function, f(r) = 1 for 
Ti < r < Tf and f(r) =0 otherwise. We further simplify 
calculations assuming the isotropic emission: X^ e \l) = 
l( e \ Then we get 



Fi°\t) 



2ttX^ 



dTf(T)l 



dl_ 
dr~ 



(9) 
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FIG. 1: The dotted lines are the emission points for the same 
arrival time parameter. 'a' and 'b' denote the points of the first 
lights with / = from Ri and Rf respectively. The solid line 
from 'c' to 'd' are the points where the last tangential ray for 
each R, Ri > R > Rf, is emitted, 'c' denotes where the light 
is emitted tangent ially from Ri. 'd' is the point from which 
the last ray is emitted from Rf. 



lating the arrival time differences, l\dl/dr\ and bending 
angles. The arrival time difference between the first ray 
and the last ray is calculated to be AT = 42. OM, while 
the corresponding proper time interval is At = 3.6 M. 
The invariant velocities are vj = —0.38 and —0.52 for Ri 
and Rf, respectively. The first ray from the surface, Rf, 
arrives at Sf = 0.17 and the last ray from Ri arrives at 
St = 0.37. One can see in Fig. [2] that the maximum of 
the flux is observed around S ~ 0.2 as discussed above. 




The integration over r for a given t is equivalent to the 
integration along the dotted line in Fig. [T]for the same ar- 
rival time difference. The duration of the observed light, 
AT, is determined by the difference between the arrival 
time of the last ray from Rf emitted at the point 'd' and 
the arrival time of the first ray from Ri emitted at 'a' 
with I 0: 



AT 



At_(l T ,rf,r l )=&\r f )-&\r z ) 
2T(Z T , r t ) - T(Z T , R(r)) +Ri + R(r f ) (10) 
(R(r f )-2M)(Ri-2M) 



2r t + 2M In - 



(r t - 2Mf 



In what follows it is convenient to use a dimensionless 
time parameter 5, which is a normalized arrival time dif- 
ference defined as S = At±/AT. The time parameter 
S changes in the interval [0,1]. For a given radius R, 
the time parameter for forwardly emitted light increases 
as I increases from I = to l ma x(R)- Then the back- 
ward emissions takes place for I < l m ax(R) an d ends 
with It(R) when 5=1. 

For the radial rays (I = 0) emitted from 'a' and 'b' in 
Fig. 1 the time parameter S are and Sf, respectively. 
Denote by St the time parameter for the point 'c' (the 
last ray from Ri, see Fig.l). Then during the interval 
[8f, St] the distant observer receives light emitted in the 
radius domain [Ri,Rf]. In other time intervals, earlier, 
< S < 5f, or later, St < S < 1, only part of this 
radius domain contributes. For this reason it is natural 
to expect the maximum flux at some S between Sf and 
St- 

To illustrate characteristic features of the light curves, 
we plot, in Fig. 2, the observable flux for the case when 
the collapse starts at the radius Ro = 9.0M and the emis- 
sion takes place from Ri = 6.0M to Rf = 4.6M. We use 
the analytic approximation developed in [10[ for calcu- 



FIG. 2: The flux from Ri = 6.0M through R f = 4.6M in 
arbitrary unit. The horizontal axis is the normalized arrival 
time S. 



Let us discuss now the time dependence of the observed 
redshift factor. For a moment, we consider the radiation 
of a sharp-in-time profile emitted at the moment of r 
when the radius is R e . Figj3]shows the redshift factor as a 
function of S for R e = B, = 6.0M, R e 5.3M, and R e 
Rf = 4.6M. One can see that for all three curves have 
the same minimum value of the redshift factor, &t = 
$(R,l T ) = 1-13, at S = 1. As discussed in Ref.[10], the 
basic reason is that the redshift factor of the last ray 
(I = It) for a freely falling surface does not depend on 
the radius of radiation but depends only on the initial 
radius Rq: 



$7 



1 



y/i - 2M/R 



(11) 



For the emission during the finite interval of time, the 
lights rays which arrive at the time S are emitted at dif- 
ferent radii, and hence have different redshift s. For each 
value of S (except for S = and S = 1) one has a finite 
range of possible redshifts and the redshift curves are 
broaden. The spectral broadening of the monochromatic 
radiation is a function of arrival time. 

The shadowed region in Fig. [H shows the broadening of 
redshift factor as a function of the time parameter S. We 
again choose Ro = 9M. The redshift factor of the first 
ray (emitted from 'a' in Fig. [Q) is <I> = <!>(/? = 6.0M, / = 
0) = 1.82. The maximum redshift occurs for the ray 
emitted from 'b' in Fig. [lj which arrives to the distant 
observer at Sf = 0.18. The maximum redshift factor 
<&max $(R 4.6,/ = 0) = 2.36 in Fig. [3] corresponds 
to the spectral broadening at Sf in Fig. [U From St = 



4 




5 

FIG. 3: The three lines 1, 2 and 3 are the redshift factors 
for the radiation emitted from R = 6.0M, 5.3M and 4.6M 
respectively. Arrival time differences are those from sharp-in- 
time pulses fTol ] . 



0.37 to S = 1, the distant observer receives the last rays, 
which are emitted tangentially for given radii, with the 
same redshift factor, $>t- Hence the spectral broadening 
for St < S < 1 has a constant minimum value of &t 
between c and d, as shown in Fig. [H 
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FIG. 4: The broadening of red shift factors. 'a,b,c' and 'd' 
are those in Fig. [T] 



In our model, the physical parameters characterizing 



the collapsing object are its mass M and three of di- 
mensionless parameters Ro/M, Ri/M and Rf/M, which 
determine the initial radius and the initial and final ra- 
diation emitting radii, respectively. In case we do not 
know the frequency of the emitted radiation v we cannot 
determine the redshift factor directly by observing the 
spectrum of the radiation. Nevertheless, if it is possi- 
ble to determine the frequency y\°^ st of the last ray with 
sufficient accuracy, then the relative (normalized) red- 
shift factors 4> = = v ^ l v \°ast can ^ e determined. 
Then by measuring Sf and St as well as <£>o and <& m ax 
one can determine Ro/M^Ri/M and Rf/M. Once we 
know Rq/M, <$>t in Eq.(fTT]) can be evaluated and we can 
determine the original frequency of emission as given by 

v = $t v \°asv ^ ne mass °f the collapsing object can 
be inferred from the observed value of At to complete 
the determination of the physical characteristics of the 
collapsing object. 

To summarize, we discuss in this work the light curves 
and the spectral broadening of the radiation emitted from 
a surface of a collapsing object. We demonstrate that the 
spectral broadening occurs when the radiation from the 
collapsing surface takes place during a finite duration of 
time. It is because of the variance of the frequency shifts 
of the light rays subjected to the gravitational redshift 
and the Doppler shift of the collapsing surface. In a sim- 
plified model of monochromatic radiation from a freely 
falling spherical surface, we discuss the possible way how 
to infer the physical parameters, the mass, the radii of 
the emission and the frequency of the radiation from the 
light curves and spectral broadenings. 
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